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SUNMARY 


This report describes the preliminary development of a N.ASTRAN model of 
the B-1 aircraft No. 2 (A/C- 2). The NASTRAN model is being developed as part 
of the airloads research study (ARS). The basic intent of the overall program 
is to utilize flight data acquired during B-1 flights, perform analyses of 
these data beyond the scope of Air Force requirements, and prepare research 
reports that will add to the technology base for future transport aircraft. 
Efforts are scheduled in distinct tasks with separate reports for each task. 


During this contract phase, a review was conducted of A/C -2 internal loads 
models to determine the minimum model complexity necessary to fulfill all of 
the airloads research study objectives. Typical model sizings were tabulated 
at selected vehicle locations, and scale layouts were prepared of the NASTRAN 
structural analysis model. The structural stiffness data, mass properties 
data, and NASTRAN model layouts are included in the reference reports, rather 
than in this document. 


INTRODUCTION 


A/C-2 (figure 1) is being employed in the airloads survey flight test 
program. This aircraft has undergone extensive ground testing to calibrate 
the strain gages utilized in the airloads survey. The aircraft provides a 
reasonable simulation of a future transport aircraft since it employs the 
large flexible structure (figure 2) envisioned in future transport designs. 
The task described in this report is part of the .ARS. 


The airloads data gathered during the flight test program can be utilized 
in the evaluation of NASA computer programs recently developed to enhance the 
analytical techniques of predicting aeroelastic response of large flexible 
aircraft. These analytical techniques include computerized structural analysis 
programs such as NASTRAN and FLEXSTAB. 
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B- i aircraft. 
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I'igiirc 2.- St ructiiral breakdown. 




Since the B-1 development program involves all experimental tests needed 
to correlate the analytical predictions with actual measured results, detailed 
plans for constructing a NASTRAN structural model of the B-1 airframe, suitable 
for use on the NASA/DFRC Cyber computer, were initiated. This model is of 
minimum complexity to give satisfactory flexibility characteristics for 
FLEXSTAB aeroelastic analysis. Included in this model are the control sur- 
faces, the control system stiffnesses, and the secondary loading edge and trail- 
ing edge structure. The plans encompass DMAP programming required to produce 
FLEXSTAB level 1.02.00 input data. 

A/C- 2 internal loads models were reviewed to scope the NASTRAN model to 
the minimum number of nodes, elements, and degrees of freedom for each struc- 
tural component. The NASTRAN model will consist of seven substructures which 
will provide capabilities of instituting structural changes on any one struc- 
tural component. A set of scaled model layouts showing the nodal point loca- 
tions for each substructure was prepared. Included in these drawings are 
typical structural sizings of the structure employed in A/C-2. 


AIRCRAFT DESCRIPTION 


The B-1 aircraft is a prototype long-range supersonic bomber with the 
capability of high-speed flight at low altitude. Configuration dimensions and 
general arrangement are presented in figure A-1. The aircraft utilizes a 
blended wing-body concept with variable -sweep wings, a single vertical stabi- 
lizer with a three-section (upper, intermediate, and lower) rudder, and hori- 
zontal stabilizers which operate independently to provide both pitch and roll 
control. The variable -sweep (15 to 67.5 degrees) wing, equipped with slats, 
spoilers (which also function as speed brakes), and flaps, provides the air- 
craft with a highly versatile operating envelope. Canted vanes, mounted on 
each side of the forward fuselage, are part of the structural mode control 
system which reduces structural bending oscillations in the vertical and 
lateral axes. 

The aircraft is powered by four YFlOl-GE-100 dual-rotor augmented turbofan 
engines in the 30 ,000 -pound- thrust class. The engines are mounted in twin 
nacelles below the wing, approximately at the left and right wing pivot points. 
For supersonic speeds, an air induction control system varies the internal 
geometry of the nacelle inlet ducts to maintain the required airflow to the 
engines for all flight conditions. 
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Fuel is carried in integral tanks in the fuselage, wing carry- through, 
and wing outer panels. The fuel system is pressurized and inerted by nitrogen. 
Fuel transfer sequencing is automatic and provides center -of- gravity control. 
The aircraft has both in-flight and single-point refueling capabilities. 


Fuselage 

The fuselage (figure A- 2) is constructed primarily of aluminum alloy 
materials arranged in a semimonocoque skin- frame -longeron type of construction. 
Titanium is used in the wing carry- through structure, in the nacelle and tail 
support structure, for various other structures where high load concentrations 
exist, and on aft fuselage skins v>/here high temperatures and acoustic levels 
are prevalent. Dielectric materials such as polyimide quartz and fiberglass 
are used for radomes and antenna covers. 

The fuselage structure is fabricated in six major sections and then mated 
together prior to attaching wings, empennage, landing gears, and nacelles. 

The following functional description of each section will provide a better 
understanding of the overall fuselage and its relationship to most of the 
aircraft subsystems. 

The crew module assembly provides a sealed enclosure with crewmember pro- 
visions and is an ejectable unit for emergency escape. The structure is 
capable of pressurization for a 2439-meter (8000 feet) altitude environment 
and incorporates a dear-vision windshield designed to bird-proof requirements, 
additional crew windows, an entry door, and an emergency exit hatch for ingress 
and egress. The floor structure supports crew seating and ejection rocket 
loads. An unpressurized section aft of the crew quarters houses the escape 
system parachutes and provides support for the stabilizing fins. Two sets of 
deployable mechanical stabilizing spoilers are hinged in the side panel frame- 
work and at the lower forward edge of the module. Structural ties to the 
forward fuselage are severed by explosive charges for emergency escape. 

The forward fuselage section includes the nose radome, a forward avionics 
compartment, an in-flight refueling receptable, the nose gear well and support 
structure, a central avionics compartment, a section of the forward fuel tank, 
a Doppler radome, an environmental control system equipment bay, and the crew 
entry stairladder structure and mechanism installation. The section also 
includes many other items of equipment such as antennas and pressure sensing 
devices. Left and right structural mode control fin surfaces are mounted on 
this section. Many large and small access doors are provided due to the high 
density of equipment installations in this assembly. 
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The forward intermediate section houses the forward and center weapons 
bays. Major bulkheads between the two bays and at each end of the bays provide 
support for the rotary weapons launchers. The aft bulkhead location also forms 
a part of the wing carry- through section. The forward bulkhead is also a 
closeout for the forward fuselage section. Large integral fuel tanks are 
incorporated into the forward intermediate fuselage structure, immediately out- 
board of the weapons bays. A systems routing tunnel occupies the upper struc- 
ture area between longerons. Provisions for avionics are incorporated in the 
side fairing area, consisting of equipment bays, antennas, and radomes. Pro- 
visions for external stores pylons, wing sweep actuation components, and flap 
and slat drive mechanisms are also incorporated in the forward intermediate 
fuselage section. 

The aft intermediate fuselage consists of the main gear well and the aft 
weapons bay. It incorporates a flight controls mixer compartment and a fuel 
tank above the main gear well. The gear uplock support structure is in the 
mixer compartment. Avionics provisions are- made in the compartment between 
the wheel wells and in the structural conpartments outboard of the wheels. 
Bulkheads at the forward and aft ends of the weapons bay support the weapons 
bay rotary weapons launcher. As in the forward intermediate fuselage, fuel is 
stored outboard of the weapons bay. A double -support frame for the aft portion 
of the nacelle extends outboard to the centerline beam of the nacelle. This 
support is approximately midway between weapons bay bulkheads. The upper 
centerline longeron and the lower outboard longerons are located and constructed 
so as to provide a high stiffness -to-weight ratio. The upper centerline 
longeron extends forward into the wing carry- through section and aft to the 
vertical stabilizer front spar. 

The aft fuselage is a semimonocoque structure and consists of the aft fuel 
tank area, the dorsal area, the aft avionics bay, and the tail cone. The tank 
area is closed in the forward and aft ends by bulkheads. The forward bulkhead 
separates the aft fuel tank from the aft weapons bay. The aft bulkhead closes 
the tank and provides mounting support structure for the horizontal tail spindle 
fitting and the aft avionics bay. The dorsal area is a dry tunnel space which 
houses flight control cables and hardware and provides for the routing of the 
electrical conduits. 


Wing 

The \>/ing consists of the wing pivot, outer wing panel, flaps, slats, and 
spoilers (figure A-3). The wing pivot consists of the pin, bearings, and 
inboard and outboard lugs, with provisions for attachment to the wing carr>'- 
through fuselage section and the wing outer panel. 
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The wing outer pcinel consists of a structural box with leading edge slats, 
trailing edge flaps, :ind spoilers over the flap leading edge. The outer wing 
is mounted on pivot bearings whose supporting lugs are mechanically attaclied 
to the wing covers. Provisions for integral fuel containment is provided in 
the outboard wing structural box. Access is provided for sealing, inspection, 
servicing, and replacement of fuel system components. Control surfaces on the 
wing include flaps, slats, and spoilers. 

The flaps are aft of the wing rear spar and are mounted on rollers between 
curved tracks. Flap actuating jackscrews are in the mid bay of the flap panels. 
Segmented leading edge slats are provided. Each segment is supported on tracks 
mounted on rollers attached to the fixed leading edge structure. Segmented 
wing spoilers are aft of the wing rear spar and above the flaps. 


Nacelle 

The nacelle is constructed of aliuninum alloy, titanium alloy, and stain- 
less steel and fiberglass laminates. Structural tvT^e is semimonocoque with 
skins, frames, longerons, and a honeycomb sandwich duct (^figure A-4^. Each 
nacelle is fabricated in two major sections, the fonvard section and the 
engine compartment. 

The forward section consists of the inlet section, the duct assembly, the 
ranps, and the center beam. The inlet section consists of the center splitter 
wedge and the upper and lower leading edges. Portions of the upper and lower 
leading edges are porous for boundary layer control. The duct assemblies 
consist of engine air intake ducts supported by frames and stringers and 
covered with an external skin. In the forward area, the duct wall is covered 
with aluminum machined skin. The intermediate and aft duct walls are cov'ered 
with fiberglass honeycomb sandwich. The inboard wall of the duct is made up 
of a fixed duct and movable ramps which provide for a variable-geometrv' system 
for air induction control. The center beam consists of four main longerons, 
interconnecting shear panels, appropriate frames, and the forward nacelle 
attach point, and is the primary vertical bending member of the nacelle. 

The engine compartment consists of the principal firewall bulkhead, the 
structure between the two engines, the primary engine attach points, and the 
aft nacelle attach points. Large hinged engine access doors are provided to 
complete the engine enclosure. Construction is frame, skin, and longeron. 

The nacelle is attached to the aircraft at four points. The fonvard 
attach point is a single fitting on the top of the centerbeam structure, which 
is connected to the wing pivot pin through a ball joint. The other three 
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attach points are in the engine section in line with the rear engine support. 
They consist of links, two vertical and one horizontal, which connect the 
nacelle to the heavy support frame extending from the aft intermediate fuselage. 


Horizontal Stabilizer 

The horizontal stabilizer (figure A-5) consists of left and right slab 
panels attached to a steel spindle projecting out of the aft fuselage stub 
structure. Both left and right panels rotate on bearings and are independently 
controlled in order for the stabilizer to provide pitch and roll control of the 
aircraft. Each panel consists of a main structural box, a leading edge assembly, 
a trailing edge assembly, a tip assembly, an aerodynamics chord plane seal at 
the inboard end, and an air seal around the spindle. 


Vertical Stabilizer 

The vertical stabilizer consists of the main box structure, leading edge 
assembly, tip assembly, and trailing edge structure (figure A-6) . The main 
box assembly supports the two upper rudder segments. Routing tunnels are pro- 
vided in enclosed areas of the main box structure for electrical and cooling 
lines required to support avionics and antenna equipment in the tip and leading 
edge components. The rudder consists of three segments. The upper two seg- 
ments are attached to the vertical stabilizer through power hinge fittings and 
actuated by hydraulic motors in the horizontal stabilizer actuator fairing. 

The lower rudder segment is supported by conventional hinge fittings and 
actuated by linear actuators between the rudder and aft fuselage structure. 

The vertical stabilizer is attached to the aft fuselage principally through 
a double-shear attachment provided on the horizontal stabilizer spindle fitting. 
The vertical stabilizer is mechanically attached to the spindle fitting by close- 
tolerance, high-strength bolts. 


NASTRAN MODEL PLANS 

Detailed plans for the finite-element modeling of the A/C-2 structure, 
intended for use with the NASA/COS'lIC release of NASTRAN level 16.0 on the 
DFRC CDC Cyber computer were prepared. 

The detailed plans constrained the model to the minimum complexity to give 
satisfactory flexibility characteristics for FLEXSTAB aeroelastic analysis. 
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Included in the model plans are the control surfaces and provisions for their 
control system stiffnesses, as vvell as the secondar>’ leading edge and trailing 
edge structure. The NASTRAN model plans specify seven substnictures consisting 
of the following: 

(,1) Wing outer planel, flaps, slats, and outboard transition lugs 

(2) Forward fuselage structure 

(3) Aft fuselage structure, wing carr)"- through structure (WCTS) , and 
inboard transition lugs 

(4) Horizontal stabilizer, leading edge, and trailing edge 

(5) Vertical stabilizer, leading edge, and rudders 

(6) Overwing and underwing fairings 

(7) Nacelle structure 

In addition to modeling the B-1 A/C-2 airframe stnicture to represent the. 
flexibility characteristics, the model was constnicted to provide stress data 
at the airload survey strain gage locations for each component. In these 
regions, the model complexity was increased to provide the desired accuracy. 

In some regions, the complexity was dictated by the NASTRAN aspect ratio con- 
straints. hfowever, after considering all of these requirements and with con- 
siderable engineering restraint, the goal of utilizing less than 5600 grid 
points appears to be achievable. 

The principal elements selected for use in the model are rods , bars , shear 
panels, isoparametric membranes, and plates. Tlie fuselage shell consists of 
rods and shear panels, with isoparametric membranes used in selected regions 
where high aspect ratio dictates their desirability. The fuselage frames and 
bulkheads will utilize bars, rods, and shear panels. The wing and empennage 
structural model will consist of bending plates, bars, rods, shear panels, 
and membranes. The use of these elements and their distribution are shown 
in table I. 

The model plans are described in the scaled model layouts showing the 
nodal point locations for each vehicle component and t>'pical stnactural 
sizings for each component. 

Detailed plans for the wing lower and upper surfaces are presented in 
figures B-1 and B-2. Wing detail views are shown in figure B-5. 
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TABLE I.- PRELIMINARY STATISTICS FOR MS NASTRAN MODEL 



Description 

o£ 

conf)onent 

NASTRAN model elements 

Substructure 

No. 

No. of 
grids 

Rod 

Shear 

panel 

Membrane 

Plate 

Bar 

1 

Wing , flaps , 
outboard lug 

1000 

550 

275 

775 

175 

400 

2 

Fwd fiiselage 

350 

450 

350 

150 

50 

300 

3 

Aft fuselage, 

WCTS , inboard lug 

650 

600 

700 

250 

275 

550 

4 

Horizontal 

stabilizer 

540 

775 

400 

450 

50 

10 

5 

Vertical 

stabilizer 

340 

325 

200 

200 

50 

10 

6 

Fairings 

300 

150 

50 

30 

200 

100 

7 

Nacelle 

340 

650 

150 

50 

250 

500 


Total 

3520 

3500 

2125 

1905 

1050 

1870 


Element 


NASTRAN nomenclature 

Rod 

= 

CONROD 

Shear panel 

= 

CSHEAR 

Membrane 


CQDNEM 1 

Plate 

= 

CQUAD 1 and CQUAD 2 

Bar 

=. 

CBAR 
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Plans for the forward fuselage substructure depicting the side, top, and 
bottom views are presented in figures B-4, B-5, and B-6. Details of the for- 
ward fuselage model are shown in figure B-7. 

Aft fuselage model substructure plans depicting the side, top, and bottom 
views are presented in figure B-8, B-9, and B-10. .Aft fuselage model details 
are shown in figure B-11. 

Horizontal stabilizer model substructure plans are shown in figure B-12. 
Detail views of the horizontal stabilizer model and actuation system are pre- 
sented in figure B-13. 

Vertical stabilizer model substructure plans, including the leading edge, 
bullet fairing, and three rudders, are presented in figure B-14, 

Plans for modeling the overwing and underwing fairing structures are 
presented in figure B-15. This complex structure includes the lower pivot 
fairing, upper pivot fairing, overwing movable fairing, forward intermediate 
fixed fairing, and underwing intermediate and aft fairing panels. 

Plans for modeling the nacelle structure, including the engine access 
doors, are depicted in the views of the top nacelle, bottom nacelle, inboard 
nacelle, outboard nacelle, nacelle centerbody, and typical nacelle frames shown 
in figures B-16, B-17, B-18, B-19, B-20, and B-21, respectively. 


-ANALYSIS STRATEGY 


General Problem Flovv 

NASTRAN, as described in ref. (1), consists of a number of subprograms or 
modules that are executed according to a sequence of macro instruct ions con- 
trolled by the executive system. A number of such sequences, called rigid 
formats, are permanently stored in the program and can be selected by means of 
control cards. Each rigid format corresponds to a particular type of structural 
analysis. The user may, in addition, devise his own sequence of module calls, 
referred to as direct matrix abstract program (DNIAP) , for problems that do not 
conform to the one of the available rigid formats. 

Figure 3 shows a simplified flow diagram for basic static analysis. Each 
block in the flow diagram represents a number of program modules. The struc- 
tural matrix assembler generates stiffness and mass matrices for the grid 
point data obtained from tabular information generated by the input file 
processor and the geometry processor. The mass matrix is used in static 
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GPt, GP2, etc, are NASTRAN functional iiK)duie names. 


I'igure 3.- NASTRAN flow chart. 


analysis tor the generation ot gravity loads and inertia loads on unsupported 
stnictures. The stiffness matrix is reduced to tlie form in which it is final Iv 
solved through the imposition of single-point and multipoint constraints, and 
through the use of matrix partitioning. Load vectors are then generatetl from 
concentrated loads at grid points, pressure loads on surfaces, gravity, 
temperature, and enforced deformation, as appropriate, and are reduced to 
final form by the application of constraints and matrix partitioning. The 
solution for independent displacements is accomplished in two steps: (1) 

decomposition of the stiffness matrix [K] into upper and lower triangular 
factors, and (2) solution for (u) for specific load vectors (P) , by means of 
successive substitution into the equations represented by the triangular 
factors of [K] . Dependent displacements are determined from the independent 
displacements by means of the equations of constraint. 


Notation System 

Many of the operations perfomied in computerized stiaictural analx'sis are 
conveniently expressed in the notation of matrix algebra. In NASTRAN, as 
described in ref. 2, matrix arrays are represented by a root s>-mbol that indi- 
cates the type of physical quantity and by one or more subscripts and super- 
scripts that act as modifiers. The root symbols used in static analysis by 
the displacement method are listed in table II. .Nearly all of the matrix 


TABLE II. - ROOT Sl'MBOLS USED IN N.ASTRAN 



1 

Vector of displacement components 

lpf 

Vector of applied load components 


Vector of forces of I'eaction 

Ivf 

Vector of enforced displacements 

[K] 

Stiffness matrix 

[M] 

Mass Matrix 

[B] 

Damping matrix 

[R] 

Matrix of constraint coefficients, as in [R] |u} = 0 

[G] 

Transformation matrix, as in ju [ = [G 1 lu 1 


• m’ m’ • n> 
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operations in static analysis are concerned with partitioning, merging, and 
transforming matrix arrays from one subset of displacement components to another. 
All components of displacement of a given type, such as all points constrained 
by single-point constraints, form a vector set that is distinguished by a sub- 
script from other sets. A given component of displacement can belong to several 
vector sets. The mutually exclusive vector sets (u ), the sum of whose members 
is the set of all physical components of displacement, are listed in table III. 


TABLE III. - VECTOR SETS 


Coordinates eliminated as independent degrees of freedom by 
multiglobal constraints 

u^ Coordinates eliminated by single-point constraints 


u 


o 


Coordinates omitted by structural matrix partitioning (i-e-, 
interior coordinates) 


u 


r 


u 



Coordinates to which determinate reactions are applied in 
static analysis 

The remaining structural coordinates used In static analysis 
(points left over) 


u ■ Extra degrees of freedom introduced in dynamic analysis to 
, describe control systems, etc 

u^ ^ + u^, the set left as analysis set for phase II 

i 

1 

Ur i u + u„, unconstrained free structural coordinates 
f a O’ 


I 

i 

i 


! 

I 

I 

i 


u U£ + u , all structural coordinates not constrained by multipoint i 
constraints ' 


Ug u^ + u^, all structural grid points including scalar points 
Up j u^ + Ug, all physical points 
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In static analysis (without substructuring) , we are concerned only with 
the grid point set (u ) and its subsets. The application of constraints and 
partitioning to the stiffness matrix involves the elimination of u^^, Ug, Uq, 
and Uj. from Ug to form a stiffness matrix referred to u^. For substructuring, 
it is necessary to perform matrix reduction as show in phase II of figure 3. 


Assembly of Structural Matrices 

The element matrix generator (BIG) and the element matrix assembler (BIA) 
generate the stiffness, mass, and damping matrices for the structural model. 

For efficiency in restart, particularly when changing from statics to d\Tiamics 
problems, the structural matrices, [K^g] , [f^gg] , [MaJ » [Baa], are assembled 
by four separate executions of EMA. §MG generates tfie various*t>-pes of struc- 
tural matrices on a selective basis. A third part of the matrix assembly matrix 
operation (SMA3) adds the contributions of the general elements to the stiffness 
matrix. 


The EMG refers to the appropriate "element" routines for calculation of 
the stiffness, mass, and damping matrices for each element. The matrices for 
each element are initially generated in an element coordinate system that is 
characteristic for each element type. The element matrices are transformed 
to the global coordinate system prior to transfer to, direct access secondary 
storage. 

The EMA assembles several columns of the structural matrices at one time. 
The number of columns assembled in one operation is limited by the space 
available in main storage. The required element matrices are transfeiTed from 
secondary storage using the direct-access read operation. The completed 
columns of the structural matrices are written on secondary storage by using 
the regular NASTRAN pack routines . 

Although the matrices generated by the stinctural matrix assembler are 
symmetric, complete columns are generated and retained for efficiency' in 
succeeding matrix operations. This is necessary because all matrix operations 
are performed one column at a time and, in dynamics applications, the matrices 
are not necessarily symmetric. Moreover, the availability of symmetric matrices 
by rows or by columns is advantageous in some of the matrix operations. 

Constraints and Partitioning 

Structural matrices are initially assembled in ternis of the set, u^, of 
all structural grid points, excluding only the set, u^, of extra points* 
used in dynamic analysis. The following paragraphs describe the subsequent 
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reduction o£ the structural matrices to the set, u^ , which is the set of 
coordinates that remains after all constraint and partitioning operations have 
been performed, and which is, therefore, the first set to be evaluated in static 
analysis. 

TTie structural matrices whose assembly is discussed in the proceeding 
paragraphs are: 

(1) t^ggJ structural stiffness matrix due to elastic structural 
elements 

(^) structural damping matrix of imaginary stiffness 

coefficients 

(3) [Bgg] the viscous danping matrix due to dajiper elements 

(4) [^gg] the structural mass matrix 

The reduction procedures for the [K^^] matrix are explained in the following 
paragraphs . 


Multipoint Constraints 

The multipoint constraint equations are initially expressed in the form, 

[f^g] (^g) - 0, 

where the coefficients are supplied by the user. Tlie user also specifies the 
degree of freedom that is made dependent by each equation of constraint, so 
that the (Ug) matrix may immediately be partitioned into t\vo subsets. 



where the set, u , is the set of dependent degrees of freedom. The matrix of 
constraint coeff^icients is similarly partitioned 

[•^g] " [Rn I Rjh] > (3) 

so that equation 1 becomes 

[Rnl (V * IRml (%) ■ 0 (41 
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[Rjj,] is a nonsingular matrix. We can, therefore, form the multipoint con- 
straint matrix. 


[Gm] = - [Rm]'^ [Rn], 
so that equation 4 may be stated as 
(4n) ■ [Gm] • 


( 5 ) 


( 6 ) 


Prior to the imposition of constraints, the structural problem may be \\rritten 
as 


[Kgg] (Ug) = (Pg), 


( 7 ) 


or, partitioning in terms of the coordinate sets, u^- and Ujjj 


- 




K 

nn 

' K 

1 nm 

1 



T 

K 

run 

' K 

1 mm 

1 


' 


( 8 ) 


Bars over symbols are used to designate arrays that are replaced in the reduc- 
tion process. 

The addition of constraints to the structure requires that the forces of 
constraint be added to the equilibrium equations. The forces of constraint are 
proportional to the corresponding coefficients in the constraint equations. 
Thus, writing the equilibrium and constraint equations together in partitioned 
form. 


K IK I G ^ 
_jin I Jim 1 jn_ 

IfT I If I _T 

"^nm I ^^mm I ^ 




G I -I I 0 
m I I 




(. 9 ) 


where [q„^J is the vector of constraint forces (um) . Straightforward elimina- 
tion of u„j and qj^ gives 


[k +K G + G^K^ + G^K g1 (u) 
[nn nmm m nm mmmmj n 



( 10 ) 
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or 


[Krail (%) = (Pn) 


( 11 ) 


where 


K =K +KG+G^K^ 


+ G^ K G 


m nun m ’ 


( 12 ) 


rin nn run m m run 


and 


P = P + G^ P . 


(13) 


n n mm 


The initial partition of Kgg and the operations indicated by equations 5, 12, 
and 13 are performed by appropriate modules of the program. The multipoint 
matrix, G^, is used in structural matrix reduction (equation 12) , load vector 
reduction (equation 13) , and data recovery (equation 6) . It is saved for these 
purposes in an auxiliary storage file. The other structural matrices, [Kgg], 
[Bgg] and [Mgg] , are transformed by formulas that are identical in form 
to equation 12. 


where (Y3) is a vector of enforced deformations, any or all of whose elements 
may be zero. The set, Uj^, is partitioned into u^ and U£ (the free or 
unconstrained set) 


Single-Point Constraints 


Single-point constraints are applied to the set, u^, in the form 



(14) 



(15) 


The stiffness matrix, Knn> is similarly partitioned 



ri6) 
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The complete structural equations, including the single-point forces o£ con 
straint, Qg, may be written in partitioned matrix form as 



Straightforward elimination gives 

[Kff] (U£) = (Pf) ■ [^fsl ^^s^ ~ (Pf) • 


The forces of constraint are recovered by means of the middle row of equa- 
tion 17; i.e. , 




(19) 


Thus all three of the distinct partitions of K^n (i.e., Kff, Kfg and Kgg) are 
needed in subsequent calculations and are placed in auxiliary' storage. For 
the other structural matrices (K^^* 

are saved. The assumption is made, implicitly, that the effects of the other 
structural matrices on the single-point forces of constraint may be ignored. 


Partitioning 


At user option, the set of free coordinates, uf, may be partitioned into 
two sets, uo and Ua, such that the Uq set is eliminated first. Thus, 



( 20 ) 


The equilibrium equations after the elimination of constraints (equation 18) 
may be written in partitioned form as 


'k 'k ' 

(u ) 

aa 1 ao 




T 1 „ 

)u I 

K K 

V o) 

L ao 1 oo_ 




( 21 ) 
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Rearrange the bottom half of equation 21 

[Koo] K) = (Pq) - (^a) 

and solve for (Uq) : 


(Uq) = [Koo]'^ W - [5^ao]^ (%)• _ 

(Note that, in practice, stiffness matrices are never inverted due to excessive 
computer running time. The practical- alternative will be explained presently.) 
Substitute for Uq into the top half of equation 21. 

[Kaa- KaoKoo'^ Kao'^l (^a) = (Pa) ' [Kao] [Koo]'^ (Po) • (--^ 

It is convenient to define the matrix 

[Go! = - [Koo]‘^ [Kao]'^. - 


so that equation 24 becomes 

[Kaa KaoGo] (^a) = (Pa) [Gol'^ (Po) , 
where advantage is taken of the symmetry of [Kool* 


Following the practice of condensation established in preceding 
paragraphs , 


[Kaa] (%) = [Pa). 

where 

[Kaa] = [Kaa] ^ [^aol [^o] 
(Pa) = (Pa) " [Go]"^ (Pq) 
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The [Gq] matrix defined in equation 25 is obtained practically from the 
solution of 


[Kool [Gol = - IKaol'^. 

where [K^q]^ is treated as a set of load vectors. Each such vector produces 
a column of [Gq] . The [Koq] matrix is first decomposed into lower and upper 
triangular factors. Once (ua) is obtained, the set of omitted coordinates, 

(Uq) , is obtained as follows. Define the set (ug) as the solution of 

[Kqo] (u°) = (Po) 1311 

Note that the triangular factors of [Kqo] obtained in connection with equa- 
tion 30 are saved for use in connection with equation 51, which cannot be 
solved until the load vector (Pq) is formed. Then, using equations 25 and 51 
in equation 22, 

(Uo) = (ug) + [Go] (ua). ' (52) 


The Guyan Reduction 

The Guyan Reduction (ref. 3) is a means for reducing the number of degrees 
of freedom used in dynamic analysis with minimum loss of accuracy. Its basis 
is that fewer grid points are needed to describe the inertia of a structure 
than are needed to describe its elasticity with comparable accuracy. If 
inertia properties are rationally redistributed to a smaller set of grid points, 
the remaining grid points can be assigned to tlie Uq set described in the 
preceding paragraphs and eliminated, leaving only the smaller u^ set for 
dynamic analysis. 

In the Guyan Reduction, the means by which inertia and damping are redis- 
tributed are to consider the [Gq] matrix of the preceding paragraphs as a set 
of rigid constraints, such that 

The [Gq] matrix now has the same implications for the Uq coordinates that 
the multipoint constraint matrix, [Gjjj] , has for the Um coordinates (equation 6). 
The reduced structural mass matrix is, by analog)^ with equation 12, 

[Maa] = [Maa ^ Mao Gq + Mao + gJ Mqo Gq] • (541 
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The reduced damping matrices, [K 3 ^] and [B^^] , are formed in the same 
manner. The structural stiffness matrix [K^^] , is given by equation 28. The 
reduced dynamic load vactor is, by analog)" with equation 13. 

(Pa) = (Pa ^ gJ P^) . (551 

The approximation made in the Guyan Reduction is that the (u°) in equa- 
tion 32 is neglected; i.e. , that the deformations of the Uq set relative to 
the Ua set due to inertial and other loads applied to the Uq set are neglected. 
The -error in the approximation is small, provided that the (u^) set is judi- 
ciously chosen. The selection should be based, in part, on an estimate of the 
relative deformation, (u^) . Thus, the members of u^ should be uniformly dis- 
persed throughout the structure and should include all large mass items. The 
basic assumption made in the Guyan Reduction is identical to that made in 
forming consistent mass matrices for individual elements. 


DMAP PROGRANMING PLANS FOR FLEXSTAB 


Detail plans for DMAP programming of the NASTRAN model to produce the 
input data required for FLECSTAB Level 1.02.00 will be accomplished during the 
structural component model generation. The development of these plans includes 
the selection of the grid points to be retained for use in the FLEXSTAB model . 
It is planned that the initial selection of points will be reviewed by NASA 
DFRC and AFFDL technical personnel for conformation. These final selected 
points will then be utilized using the CMIT and/or ASET cards to achieve the 
matrix reduction, with the suitable matrix being generated within the NASTRAN 
computer program processing. The desired flexibility matrix can be generated 
for use by FLEXSTAB by using the OUTPUT2 module in conjunction with the SOL\T 
module of NASTRAN. Since there has been limited application of N.ASTRAN 
Level 16.0 in conjunction with FLEXSTAB Level 1.02.00, the technical personnel 
from Rockwell, NASA DFRC, and AFFDL will jointly develop the detailed plan 
for the DMAP programming. 
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Figure A-1. - General arrangement - RDT§E A/C-1 and-2 
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GENERAL ARRANGEMENT 
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Figure A-1. - Continued. 
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1. PITOT STATIC BOOM (WITH AOA 

AND SIDESLIP VANES) 

2. FORWARD RADOME 

3. AERIAL REFUEL RECEPTACLE 

4. PITOT STATIC PROBE ♦ 

5. TOTAL TEMPERATURE PROBE ♦ 

6. STRUCTURAL MODE CONTROL SYSTEM 
VANE ♦ 

7. angle OF ATTACK VANE ♦ 

8. CREW ENTRY WAY 

9. FORWARD FUSELAGE FUEL TANK 
(TANK NO. 1) 

10. FORWARD INTERMEDIATE FUSELAGE FUEL 
TANK (TANK NO. 2) 


35. LOWER RUDDER 

36. AFT RADOME 

37. AFT AVIONICS COMPARTMENT 

38. LN2 DEWAR 

39. ENGINE NACELLE ♦ 

40. MAIN LANDING GEAR ♦ 

41. AERIAL REFUEL/WING INSPECTION LIGHT* 

42. ENTRY DOOR 

43. ENTRY LADDER 

44. NOSE LANDING GEAR 

45. LANDING/TAXI LIGHT 

46. LANDING LIGHTS (2) 


47. FORWARD AVIONICS COMPARTMENT 

48. EJECTABLE CREW MODULE 

49. FORWARD CREW STATIONS 

50. CREW SEAT (4) t 

51. ESCAPE HATCH (SEVERABLE) 

52. AFT CREW STATIONS 

53. CONTROLS FOR ENTRY LADDER. APU, AND 
MAIN GEAR DOORS 

54. SURVIVAL EOUIPMENT 

55. SIDE WINDOW (SEVERABLE) ♦ 


11. MAIN FUEL TANKS 

12. MAIN WHEEL WELL EOUIPMENT (INTERMEDIATE 
AVIONICS) COMPARTMENT 

13. AFT INTERMEDIATE FUSELAGE FUEL TANK 
(TANK NO. 3) 

14. AFT FUSELAGE FUEL TANK (TANK NO. 4) 

15. HORIZONTAL STABILIZER ACTUATOR ♦ 

16. HORIZONTAL STABILIZER 

17. FLIGHT CONTROLS MIXER BAY 

18. ENGINES ♦ 

19. FLAPS (6) ♦ 

20. SPQILERS/SPEED BRAKES (4) ♦ 

21. FUEL JETTISON OUTLET ♦ 

22. POSITION LIGHT ♦ 

23. SLATS (7) ♦ 

24. WING FUEL TANK ♦ 

25. APU* 

26. HYDRAULIC RESERVOIRS* 

27. INLET RAMP MECHANISM* 

28. WING PIVOT 

29. SUPPLEMENTAL POSITION AND 

anticollision light * 

30. WING GLOVE AVIONICS COMPARTMENT * 

31. CENTRAL AVIONICS COMPARTMENT 

32. VERTICAL STABILIZER 

33. TAIL/ANTICOLLISION LIGHT 

34. UPPER AND INTERMEDIATE RUDDERS 


♦ Both Sides (L and R) 
f Right aft seat temporarily 
removed 
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Figure A-3. - Wing- structure diagram, outer (-55B). 
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Figure A-6. - Vertical stabilizer structural arrangement. 
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Figure B-1. - Wing lower surface. 
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Figure B-2. - Wing upper surface. 
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Figure B-4. - Forward fuselage" side view. 
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Figure B-5. - Forward fuselage top view. 
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Figure B-6. - Forward fuselage bottom view. 
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Figure B-7. - Forward fuselage details. 
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Figure B-8. - Aft fuselage side view. 
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Figure B-9. - Aft fuselage top view. 
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Figure B-10. - Aft fuselage bottom view. 
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Figure B-12. - Horizontal stabilizer. 
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Figure B-13. - Horizontal stabilizer detail views. 
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Figure B-14. - Vertical stabilizer. 
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Figure B-15. - Wing fairings. 
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Figure B-20. - Nacelle centerbody. 
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Figure B-21. - Typical nacelle frames. 
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